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> : ABSTRACT 
in 

We investigate a possibility of observing CP asymmetries in the partial widths for 
the decays B~ — ► 7r + 7r~ K~ and B~ — > K + K~K~ produced by the interference of the 



non-resonant decay amplitude with the resonant amplitudes. The resonant states which 



subsequently decay into 7r + 7r~ and K + K~ or K~n + are charmonium cc states with J p = 
+ , 1~, 1 + or the 0(1020) meson. We find that the largest partial width asymmetry comes 
from the Xco resonance, while the resonance ip(2S) gives a partial width asymmetry of 



the order 10%. 



1 INTRODUCTION 



The experimental data on B mesons decays into three mesons accumulate jl] - 0j and a 
number of important questions on their decay dynamics and their relevance for the precise 
determination of the CP violating phase 7 should be answered |Sj - [Hj. Motivated by 
Belle and BaBar results on the B mesons three-body decays jU El El Ej we continue 
with the study of CP violating partial width asymmetry in the B ± — > K ± 7i + n~ and 
B^ — > K ± K + K~ decay amplitudes. 

Recently, we have studied a case of the partial width asymmetry resulting from 
the interference of the non-resonant B~ — > M + M~K~ , M = ir,K, and the resonant 
B~ — > XcqK~ — > M + M~K~ decay amplitudes In both decay modes, the dominant 
contribution to the non-resonant amplitude comes from the penguin operators. However, 
there is a small tree level contribution in which enters the weak CP violating phase 7. 
The strong phase, which is necessary to obtain the CP violating asymmetry, enters trough 
the dispersive part of both non-resonant and resonant amplitudes. 

It was pointed out by the authors of and (THj that the dispersive part of the non- 
resonant amplitude exactly cancels the dispersive part of the resonant amplitude coming 
from the intermediate state which is identical to the final state. Therefore, the partial 
width asymmetry for i? ± — > 1ZK ± — > M + M~K ± , M = tt, K, will be proportional to 
the decay width of the resonant state 1Z to all channels excluding the M + M~ state. It 
means that one would expect a large CP asymmetry for the two-meson invariant mass 
in the Xco mass region since the decay width of Xco is rather large and its branching 
ratio to M + M~, M = ir, K is negligible. The amplitude for the Xco resonant decay 
mode was determined using the narrow width approximation jSJ [7] and the experimental 
results for the B~ — > XcoK~ and Xco M + M~ decay rates. The asymmetry was 
found to be about 20%. In the case of B~ — > K~ M + M~ there are, however, additional 
important reasons why the partial width asymmetry can be sizable. In fact, if in the 
B~ — > K~ M + M~ decays the partial widths coming from the non-resonant Ai m and the 
resonant Ai r amplitude are of the same order of magnitude, as in our analysis at the Xco 
resonance region Jo], one obtains a significant CP violating asymmetry. In the case of 
negligible non-resonant amplitude relative to the resonant amplitude (or vice versa) one 
would get a very small partial width asymmetry. 

In this paper, we extend this analysis to the case of the CP violating partial width 
asymmetry when the interference with the non-resonant amplitude occurs in the neigh- 
borhood of the resonance 7Z which is a charmonium cc state with J p = + , 1~, 1 + or a 
light vector and scalar meson. We will restrict our investigation only to those resonant 
states 1Z for which the decay B~ — > TIM", M = K, 7r amplitude does not have two or 
more contributions with different weak phase, as from the experimental branching ratio 
we are able to extract only the absolute value of the amplitude. For example in the case 
of B~ — > 1ZK~ with 1Z = p° there is a penguin and a tree amplitude and one needs to 
know their relative sizes to constrain the partial width asymmetry. In this decay mode it 
has also been found that the naive factorization fails to describe the decay rate [TBI [T7] . 
Therefore, we concentrate on the partial width asymmetry for the cases in which the 
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relevant two-body amplitude can be completely extracted from the measured decay rates. 

In the case of the B~ — > if ~7r + 7r~ partial width asymmetry, the intermediate reso- 
nant states of interest would be the light strange mesons if* (890), if i (1270), if i (1400), 
ifo(1430) etc. in the decay chain B~ — > TZtt~ — > K^-k + -k~ and the charmonium cc states 
in the decay chain B — > 1ZK~ — > if~7r + 7r~. The B~ decays to these strange mesons in 
the final state occur as a pure penguin transition. Among all such decays only the rates 
for TZ = if* (890) and if* (1430) were measured |3]. However, the if* (890) and if* (1430) 
mesons decay to if ~7r + with the branching ratios close to 100%. In the case which we 
consider it means that the partial decay width to the rest of the states is negligible and 
the corresponding CP violating asymmetry vanishes. The relevant charmonium cc states 
in the decay chain B~ — > 1ZK~ — > tt + tt~K~ are produced by the b — > ccs transition. 
The resonant B~ — > M + M~K~ amplitude is obtained from the tree level contribution 
which is proportional to the V c b and V cs CKM matrix elements, followed by the strong 
decay of the cc state into 7r + 7r~ or K + K~ via the OZI (Okubo-Zweig-Iizuka) suppressed 
strong interaction. Apart from the already mentioned Xco state, this category includes 
also J/ip, Xcii Xc2, ^(25) etc. We will consider contributions from all the above mentioned 
states, even though the B~ — > Xc2-^f~ and Xci M~ M + branching ratios have not been 
measured yet. Nevertheless, we expect that the partial width asymmetry in this decay 
modes can be rather large. Although one would expect that the b — > ccs transition will 
give larger rates for the two-body decays than in the case of the b — > uus transition, the 
fact that the strong transition of the charmonium states is OZI suppressed makes the 
non-resonant and resonant partial width to be of the same size and this leads to a sizable 
CP violating asymmetry. 

In the case of the B~ —>■ K~K + K~ decays with the two-meson invariant mass below 
the charmonium production threshold, the resonant contribution comes from the inter- 
mediate ss states. We consider only the CP asymmetry at the 0(1020) resonance and do 
not consider contributions from the scalar meson resonances due to the lack of knowledge 
on their structure. 

In the analysis of the partial width CP asymmetry, one needs a knowledge of the 
non-resonant amplitudes. We compute the non-resonant decay amplitudes by using a 
model which combines the heavy quark effective theory and chiral Lagrangian, previously 
developed in jH] - [HI- This model assumes the naive factorization for the weak vertices. 
The fact that the factorization works reasonably well in the relevant two-body decay 
modes encourages us to apply it in the three-body decays we consider here. Even more, 
the experimental investigation of the non-resonant amplitudes done by Belle collaboration 
jH] indicates that one has to rely on a model when discussing the non-resonant background. 
In comparison with our previous investigation E] , we include now the contributions of 
B*(0 + ) resonances. 

In Section 2 we present the calculation and the results on the non-resonant B~ — > 
l\ M ' M . M = 7r, if decay modes, while in Section 3 we analyze the partial width 
asymmetries. The summary of our results is given in Section 4. 
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2 NON-RESONANT AMPLITUDES 



The effective weak Hamiltonian relevant for the B ± — > K ± M + M~ decays and their CP 
conjugates after Fierz reordering of the quark fields and neglecting the contribution of 
the color octet operators is [TH] - [2*T] : 

G 10 



H eff = -4(y„*K 6 (aiOi + a 2 2 ) + ^ c6 (a lc lc + a 2c 2c ) - V^^aA) , (f 
The effective Wilson coefficients are denoted by a, and the operators read: 



Oi = 2 = (u^Y-A^V-A, (2) 

Oic = (cfc)y_A(sc) v _A , 2c = (ce)v-A(sb) v „ A , (3) 
°3= Y {qq)v~A{sb) V -A, C 4 = £ (g6)v_,i(sg)v_A, (4) 

q=u,d,s q=u,d,s 

5 = E feW^y.^, 6 = -2 E (?(l-7s)6)(s(l+75)9), (5) 

q=u,d,s q=u,d,s 

3 3 

q=u,d,s q=u,d,s 

3 3 

q=u,d,s q=u,d,s 

where {q\q2)v±A stands for gi7 M (l ± 75)^. Here Oi and 2 are the tree level operators, 
O3 — 06 are the QCD penguin operators and O7 — O10 are the electromagnetic penguin 
operators. From [21] we take ai = 1.05, a 2 = 0.07, 04 = —0.043 — 0.016* and a 6 = 
—0.054 — O.OI62. The values of the other Wilson coefficients are at least one order of 
magnitude smaller and therefore we can safely neglect them. 

For the CKM matrix elements the Wolfenstein parametrization is used (V u b = AX 3 (p — 
ifj), V us = A, V u = -AX 2 , V tb = 1), with A = 0.8, A = 0.228, p = 0.118 - 0.273 (the 
average value 0.222) and fj = 0.305 — 0.393 (the average value 0.339) [22] • The matrix 
elements of the four quark operators acting in Oi for the B~ — > K~tt + tt~ decay can be 
written using the factorization assumption as: 

< 7i + 7i-K-\(sb)(qq)\B- >=< K~\(sb)\B- >< 7i~7i + \(qq)\0 > , (8) 

< 7r + 7r~K'\(db)(sd)\B~ >= 
<n-\(db)\B- >< K-7i + \(sd)\0> , (9) 

< Tt + Tt~K-\{ub){su)\B- >=< n + n-\(ub)\B- >< K~\(su)\0 > (10) 

+ < 0\(ub)\B- >< ^7t + 7t-|(sm)|0 > . 

In the above equations (qiqj) denotes the vector or axial-vector current or scalar or pseu- 
doscalar density. By analyzing the matrix elements given above, one finds |5 that only 
the first term in (fTU|) gives important contribution to the non-resonant decay rate. Terms 
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Figure 1: Feynman diagrams contributing to the non- resonant part of the amplitude. 



(jHJ and (jOJ) contribute to the resonant part of the amplitude (through resonances which 
decay into 7r + 7r~ or K~n~ respectively), while the annihilation term in (|Hlj) is found to 
be negligible as explained in jSJ. In the matrix element of the Oq operator, additional 
terms might arise, but they are either small or cancel among themselves [3]. 

The B~ — > K~K + K~ amplitude can be factorized in the same way by replacing 7r =l= 
with K ± in (jHj)- (jl(J|) . However, in this case, the contribution coming from B~ — ► p°K~ — > 
K~K + K~ (Eq. (jHJ)) is part of the non-resonant amplitude, since the p° mass is below 
the K~K + threshold. Nevertheless, we find this contribution to be small due to the 
suppression of the p° propagator in the high energy regions and due to the smallness of 
its Wilson coefficients (02 and 09) and will therefore neglect it. The same argument holds 
if the p meson is replaced by similar resonances (a etc.). 

Next, we proceed with the determination of A n = (n + (p2)ir^(pi)K~(p 3 )\Oi\B^) and 
Ak = {K + {p2)K~ (pi)K~ (p 3 )\Oi\B~) . The approach used in the calculation of these 
matrix elements was already explained in 12H1 12U • Here, we follow the same method, 
but add the contributions of the Bq scalar meson resonances. We introduce the bq states 
(q = u, d, s), with the J p = 1 + , + assignment incorporated in the field S 26 : 

S = ^(l + v al a )[D^ l5 -D ], (11) 

which then interacts with the bq J p = l - ,0~ multiplet (H) and the pseudo Goldstone 
mesons by the means of the Lagrangian: 

C s = ihTr{S bltll5 AlH a ), (12) 

where A^ 1 = 1/2 (£t<9 M £ — £<9 M £t) with the the light pseudoscalars fields in £. The weak 
current is given by 

^ = £rr[ 7 ^(l- 75 )^L]. (13) 
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The parameter h = —0.6 ± 0.2 is taken from the recent study of the D s (0 + ) state [27] . 
while for the scalar meson decay constant we use F + = 0.46 GeV 3 / 2 j2*Bj . 
The matrix element A n can be written as: 



X = -f 3 [m 2 r nr + l/2(m 2 B - m 2 - s)w™ + l/2(s + 2* - m B - 2m 2 - m 2 3 )w™)\ , (14) 

where the form factors w+, id™ and r nr are determined by calculating contributions coming 
form the Feynman diagrams in Fig. 1: 



j. 3/2 1/2 / ■> > ,\ <■ 

9 Jb*™>b* m B L m B -m{-t\ / 



+ ~ hh t-m\, [ 2m%. ) + 2 fl f 2 

v 7 ^* 2 ( 2t + 8 - m 2 - m 2 - 2m 2 ) + F+h ^ m * Z 1 
2hhm 2 B i2t + S mB ™ 3 2mi)+ 2 fl m% t - m 2 ' 



g f B mtmf f + m 2 B - m 2 - t \ + ym^ + F+fr^ m| - t 
A/a t-m\, \ 2m 2 B , I fj 2 2f x m% t - m|» 



r nr = - 9 . . /? 2T(2t + s - m| - m 2 - 2m?) + (17) 

2/1/2 (m| -m%) 2/i/ 2 

4g 2 f B m B m B * (s—m\ — m\ (t + m\ — m\){m 2 B — m\ — t)\ 



fMm 2 -m 2 B )(t-m 2 B ,)\ 2 Am* 2 
F + h^/mE m 2 B - t F + h 2 ^/m^ (m 2 B -*)(*- m|) 



+ 



j x m\ t - m 2 B * hf 2 m% (t - m 2 B0 ) (m| - m|* ) ' 

We used the Mandelstam's variables s = (p B — P3) 2 and t = (p B — Pi) 2 ■ Indies 1, 2 and 
3 correspond to 7r~, 7r + and K~ respectively (f\ — / 2 = /r, A = fx, m i = m 2 = m n , 
^3 = m-ivr)- The masses tub, mB* and mgj correspond to the £?~, B°*(l~) and S°*(0 + ) 
mesons, (1~) denoting vector and (0 + ) scalar states. The rest of parameters are taken to 
be f n = 0.132 GeV, f K = 0.16 GeV, f B = 0.175 GeV, f Bs = 1.16/ B , «i = 0.16 GeV 1 / 2 , 
a 2 = 0.15 GeV 1 / 2 as in 0. For the strong coupling g we use g = 0.56 according to the 
measurement of J2B]- Note that in ^ there are misprints in Eq.(16): the sign in front of 
a 2 is reversed, as well as the overall sign in (22). 

The matrix element of O4 has the same structure as the matrix element of 0\ while 
for determining the matrix element of 0$ we follow the approach described in Using 
the expressions (18)-(20) of [5], we find that its contribution is proportional to the matrix 
element of 0\ or O4 with the proportionality factor k§ = —2 m ^ . 

The matrix element (K + (p2)K~(pi)K~~(p 3 )\Oi\B~) is calculated in the same way. The 
expression for Ak and its form factors can be derived from the Eqs. ([14 )1 - ([17 )1 . adding 
the additional contribution obtained by interchanging s and t and by taking fi — / 2 = 
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/ 3 = fx, mi — m 2 — m 3 — n^K- I n the propagators the B meson masses are replaced by 
the B s mass. 

Now, the non-resonant amplitudes for B~ — > M + M~ K~ can be written as 

G 

M ni = -j^A K AVubV* s a x - V tb V t * s (a 4 + k 6 a 6 )) , (18) 

with An, K defined in Eq. (jl4H17jl . while for B + -> M + M-R + we have: 

M nr = -^AxAKhVusai - V; b V ts (a± + k 6 a 6 )) . (19) 

Using the above expressions, we obtain the following branching rations: 

BR(fi ± -> ^ ± 7r + 7r") nr = 9.0 x 10~ 6 , BR(B ± -> f^fT+fr ) nr = 14 x 1(T 6 , (20) 

where BR(S ± -> K ± M + M~) nr stands for the CP-averaged rates for £~ -> M~M + K~ 
and S+ -> M~M+K + ((BR(S" -> /\ .1/ .1/ ) + BR(5+ -> if + M + M-))/2). In it 
was found that due to the imaginary part of the 0,4 and ci6 Wilson coefficients we can have a 
large CP asymmetry between the non-resonant B + — > M~ M + K + and B~ — ► M~ M + K~ 
amplitudes. The size of this asymmetry depends on the p and fj CKM parameters and is 
rather large (60% in the case of B^ —>■ 7t~tt + K ± and 40% in a case of B ± — > K~K + K ± 
decay mode). The largest error in BR(B ± — > i^ ± M + M _ ) nr , due to the model parameters, 
comes from the uncertainty in the CKM weak phase 7, the decay constants and the 
coupling g. For example, by taking two times smaller g, the rate BR(5 ± — > K ± n + n^) m . 
decreases by 40% and BR(B ± -> K^+ir-)^ by 30%. Varying p between 0.118 and 0.273 
and fj between 0.305 and 0.393 gives BR( J B ± -> K^ir+ir-)^ = (6.2 - 12.6) x 10~ 6 and 
BR(5 ± — > K ± K + K~) nr = (11 — 17) x 10~ 6 . The uncertainty in the branching ratios 
coming from the B decay constants is not larger than 10%. 

The Dalitz plots for B~ — > K~ M + M~ (M = tt, K) decays, are given in Fig. 2 
(g = 0.56). We can see, that the non-resonant B~ — > K~K + K~ decay amplitude is 
rather flat, while in the case of B~ — ► 7i^tt + K~, an increase at low K and n momenta 
phase space region is evident. The inclusion of the scalar states B*(0 + ) is not giving 
significant contribution to the decay rate, increasing it by few percent in both decay 
modes. 

Recently B factories El E] got some insight into the nonresonant contribution to 
the B~ — ► K~M + M~ decay widths. The preliminary results of the Belle collaboration 
are [BE]: BR(B ± -> K+ir+ir-)^^ = 14 ± 6 x 10~ 6 and BR(B ± -»• K^+K')^^ = 
22.5 ±4.9 x 10~ 6 , while the BaBar collaboration still has only the upper limit [HH]. The 
inclusion of the nonresonant contribution in the B^ — > K 7r + ir~ Dalitz plot analysis [3 j 
was motivated by the obvious deficit of the data in the low K + tt~ invariant mass phase 
space region (see Fig. 11, first row of [3J. They used rather simple fit (see Eq. (11) 
|3j) for the nonresonant amplitude. Nevertheless, as pointed out by J. R. Fry pQ, this 
contribution is not yet well understood and more studies of this problem are expected. 
Calculated ranges for the branching ratios within our model BR(B ± —>■ K ± 7c + TT~) nT = 
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(6.2 - 12.6) x 1(T 6 and BR{B ± -> K^K+K'^ = (11 - 17) x 1(T 6 agree with the Belle 
collaboration's results within one standard deviation. Unfortunately, the experimental 
statistics is still to low to compare the distributions of the differential decay rate of 
the model and the experiment. It is interesting that our model predicts rather small 
differential decay width distribution in the region of the low tc + K~ invariant mass. In 
order to describe data given in Fig. 11 of [3] it seems that one needs such behavior of the 
nonresonant amplitude. 

In addition the results of [3] indicate existence of the broad structures in the exper- 
imental data at y/s ~ 1.3 GeV in the K + tt + tt~ final state and at y/s ~ 1.5 GeV in the 
K + K + K~ final state. Although one explanation is that light scalar resonances might be 
responsible for this effect we suggest that these increases might be induced by the 
nonresonant effects also, what can be seen in the presented Dalitz plots (Fig. 2). 



3 PARTIAL WIDTH ASYMMETRY 

For the resonances in the s-channel, the partial decay width T p for B~ —>■ MMK~ , M = 
7r + , K + , which contains both the non- resonant and resonant contributions, is obtained 
by integrating the amplitude from s min = (m R - 2T R ) 2 to s max = (m R + 2T R ) 2 : 

r p = j- ,„ QO 3 / ds dt\M nr + M r \ 2 . (21) 



(2tt) 3 32m| „ 

Similarly, the partial decay width Tp for B + — > MMK + , M = tc + , K + is defined in a 
same way. The CP violating asymmetry is then: 

A P = &^}. (22) 
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R 




TZ — > 7T + 7T 


n -> k+jst- 


J/V 

XcO 
Xcl 

^(25) 


(7.9 ±2.0) x 10" 6 
(1.01 ±0.05) x 10" 3 
(6.5 ± 1.1) x 10" 4 
(6.0 ± 2.4) x 10" 4 
(6.6 ±0.6) x 10" 4 


(1.47 ±0.23) x 10~ 4 
(5.0 ±0.7) x 10- 3 
< 2.1 x lO- 3 
(8 ±5) x 10- 5 


(49.2 ± 0.7)% 
(2.37 ±0.0.31) x 10" 4 
(5.9 ±0.9) x 10-3 

< 2.1 x 10-3 
(1.0 ±0.7) x 10" 4 



Table 1: The decay B — > 7ZK width and the branching ratios for TZ — > M + M . 

It is important to notice that at the phase space region where the invariant mass of 
M + M~ approaches the mass of the TZ resonant state, M + M~ can re-scatter trough that 
resonance as it is visualized in Fig. 3 (left figure). If Br(TZ —>■ M + M~) is large, this can 
lead to a significant absorptive amplitude and it contributes to the partial decay width 
asymmetry. As mentioned in Introduction and explained in Appendix, such contribution 
is exactly canceled by the absorptive part of a resonant decay, where the resonance re- 
scatters trough the intermediate states equal to final states (Fig. 3 (right figure)). This 
implies that one has include the factor (1 - Br (72. — > M + M~)) in the equation for the 
partial decay asymmetry. 

In the calculation of the Y p — Tp, by taking V ub = \Vub\e 11 ~ AX 3 (p — iff), we derive: 

r r . Am R T R (l-Br(TZ^M + M-)) 
T > - ^ = Sm7 (2n)*32m% 

ds I dt —=\V ub \\V* s \ai < KTnr\Oi\B > m 



x\M(B- -> RK~)\- 1 \M(R - tt-tt- 

[s - m R y + {m R v R y 



(23) 



while the T p + Tp is given by: 

r__o L_ 

{2n) 3 32m 3 B , 

Smin ^ train (^) 

x{|-MJ 2 + IwMfS" -> P-FT) =-A*0R -> tt-tt- 

s - m|j + im R T R 



(24) 



The -B — > TZK — > K M + M amplitudes are obtained from the experimental data 
and the measured branching ratios for B~ — > 7ZK~ and R — *> M + M~ are given in 
Table 1. 

For the scalar resonance exchange (xo i n our case) in the B~ — > SK~ — > M + M~K~ 
decay, we have: 

M(B~ - SP 1 {q 1 ) -> P 1 (g 1 )P 2 (g 2 )P 3 (g 3 )) = 
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K v 


9VTTTT 


9vkk{9vkit) 




2.26 x 10" 9 




6.34 


J/1P 


1.41 x 1(T 7 


2.76 x 10" 4 


4.37 x 10- 3 




2.04 x 10~ 7 


1.41 x 10- 3 


0.166 


Xcl 


1.68 x 1(T 7 


< 0.0126 


< 0.0134 



Table 2: The parameters used in our numerical calculations. 



M(B- - SP 1 (q 1 ))— J-— M(S - P 2 (g 2 )P 3 (g 3 )) , (25) 

m 23 ~ m S + 11 S m S 

where m 2 3 = (g 2 + q 3 ) 2 while mj and T$ are the mass and the decay width of the scalar 
resonance respectively. We find M{B~ — > \cqK~) = 3.34x 10~ 7 GeV, M(xco — > vr + 7r~) = 
0.118 GeV and A4(x co -> K + K~) = 0.132 GeV. 

The amplitude for the -B~ decay into light vector and pseudoscalar resonance and the 
amplitude for the vector meson decay into two pseudoscalar states are given by: 

M(B- - V(s)P 1 (q 1 )) = Kq x -e\ M(V ^ P 2 (g 2 )P 3 (g 3 )) = ^^-(<h-q 3 )-e. (26) 

The amplitude for the three-body resonant decay for this case is: 

M(B~ ^ V(e)P 1 {q 1 ) ^ P 1 (g 1 )P 2 (g 2 )P 3 (g 3 )) = 

Kg VP2 P 3 -qx ■ (gg - gg) + (gi • (g2 + qs){ml - mp)/™ 2 , 
y/2 - my + iT v m v 

where m| 3 = (g 2 + q 3 ) 2 , while m 1; m 2 , m 3 and m v are the masses of particles Pi, P 2 , P 3 
and V respectively and Ty is the width of the vector resonance. Using above formulas, 
we find the expression for the resonance exchange in the s— channel: 

M(B- - VK~ - K-M+M-) = KygvMMml + 2 m l +m ] x -2t-s 

2^2 s -ml + iT v m v 

where M stands for K or 71". In the case of the K~K + K ± mode the contributions coming 
from the s and t channels are completely symmetric. Values of Ky and gvp 1 p 2 are given 
in Table 2. 

The results for the asymmetries are presented in Tables 3-6. Tables 3 and 5 contain 
the asymmetries for g = 0.56. The off-shell mass effects might reduce this coupling as 
mentioned in and therefore we present the partial width asymmetries for g = 0.27 
(Tables 4 and 6). We calculate asymmetries for the ranges p = 0.118 — 0.273 (the average 
value 0.222) and fj = 0.305 — 0.393 (the average value 0.339) as in j22j. The subtraction 
of Br(TZ — > M + M~) in Eq. ()23j) makes a sizable effect in the case of the B~ — ► K~<p — > 
M + M~ K~ asymmetry, but it is negligible in the case of partial width asymmetry in the 
neighborhood of charmonium resonances. Then we can draw the following conclusions: In 
the case of B~ — > K~7t + tt~, all partial width asymmetries are not very large. The largest 
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p = 0.118 


p = 0.118 


p = 0.273 


p = 0.273 


p = 0.222 




fj = 0.305 


77 = 0.393 


fj = 0.305 


fj = 0.393 


fj = 0.339 




10.2% 


13.0% 


10.3% 


13.1% 


11.3% 




0.8% 


1.1% 


0.8% 


1.1% 


0.9% 




3.5% 


4.5% 


3.5% 


4.5% 


3.9% 




17.3% 


21.8% 


17.6% 


22.1% 


19.3% 



Table 3: The partial width asymmetry for B — > K tt + tv , calculated with g = 0.56 and 
given p and f/. A p (x c i) is obtained by taking the upper bound for gvMM- 





p = 0.118 


p = 0.118 


p = 0.273 


p = 0.273 


p = 0.222 




?7 = 0.305 


?7 = 0.393 


f] = 0.305 


fj = 0.393 


f\ = 0.339 


A P ^(2S)) 


13.5% 


17.3% 


13.7% 


17.3% 


15.1% 


A p (J/tP) 


1.2% 


1.6% 


1.2% 


1.6% 


1.4% 




5.0% 


6.4% 


5.0% 


6.5% 


5.6% 


Ap(Xco) 


12.8% 


16.1% 


12.9% 


16.3% 


14.2% 



Table 4: The partial width asymmetry for B — > if 7r + 7r , calculated with g = 0.27 and 
given p and f/. A p (xci) is obtained by taking the upper bound for gvMM- 





p = 0.118 


p = 0.118 


p = 0.273 


p = 0.273 


p = 0.222 




= 0.305 


?7 = 0.393 


f) = 0.305 


f] = 0.393 


f\ = 0.339 


A P (<j>) 


0.3% 


0.3% 


0.3% 


0.3% 


0.3% 


A p (i>(2S)) 


3.1% 


3.8% 


3.0% 


3.7% 


3.3% 


A p (J/4>) 


0.03% 


0.04% 


0.03% 


0.04% 


0.03% 


Ap(Xd) 


0.5% 


0.7% 


0.5% 


0.3% 


0.6% 


Ap(Xco) 


28.8% 


35% 


27.6% 


33.8% 


30.6% 



Table 5: The partial width asymmetry for B — > K K + K , calculated with g = 0.56 
and given p and fj. A p (xci) is obtained by taking the upper bound for gyuu- 





p = 0.118 


p = 0.118 


p = 0.273 


p = 0.273 


p = 0.222 




f] = 0.305 


fj = 0.393 


f] = 0.305 


f) = 0.393 


p = 0.339 


A P (<f>) 


0.3% 


0.3% 


0.3% 


0.3% 


0.3% 


A P ms)) 


8.1% 


10.1% 


7.9% 


9.8% 


8.8% 


A p (J/4>) 


0.55% 


0.71% 


0.55% 


0.71% 


0.61% 


Ap(Xd) 


3.0% 


3.8% 


3.0% 


3.8% 


3.3% 


A P (Xco) 


23.1% 


28.7% 


22.5% 


28.0% 


25% 



Table 6: The partial width asymmetry for B — > K K + K , calculated with g = 0.27 
and given p and fj. A p (xci) is obtained by taking the upper bound for gvMM- 
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asymmetry was found in the case of Xco resonance and then in the case of ip(2S). The 
partial width asymmetry A p (xci) is calculated by taking the upper bounds for x cl M + M~ 
coupling. All these asymmetries are rather stable on the variations of g. In the case of 
B~ — ► K~K + K~ the situation is different. Calculated partial width asymmetries except 
the A p (xco) are smaller than in the case of B~ — > K^-k + -k^ . They depend more on 
the variations of the g coupling. The only relatively sizable partial width asymmetry in 
addition to A p (xco) is j 4p('0(2s)). 

We have also estimated the partial width asymmetry for the B~ — > Xc2^~ channel, 
by assuming the B Xc2 k coupling to be of the same size as for the vector (scalar) mesons 
and we found it negligible. 

4 SUMMARY 

In this paper we have investigated the partial width asymmetry for the B~ — > M + M~ K~ , 
M = ir,K decays which results from the interference of non-resonant and resonant am- 
plitudes. 

First, we have calculated the non- resonant branching ratios and found that the model 
we use gives the decay rates in the reasonable agreement with the Belle collaboration 
results [3]. Comparing the Dalitz plots for the non-resonant decay modes obtained from 
our model with the experimental data we find that our model reproduces the data 
quite well. The inclusion of the Bq scalar meson is rather insignificant, contributing only 
by few percents to the rate. 

We then consider the partial width asymmetries for a few resonant decay modes for 
which the amplitude does not contain the weak phase 7. In the case of B~ — > -k + -k^K^ the 
largest partial width asymmetry arises from the interference of the non-resonant amplitude 
with the resonant amplitude coming from the Xco and ip(2S) states. In the case of B~ — > 
K~K + K~ the largest partial width asymmetry comes from the Xco scalar resonance, while 
and is about 10% in the case of ip(2S) state. 
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Figure 3: Diagrams presenting the non-resonant (left) and the resonant (right) contri- 
butions to dispersive part of the amplitude in the phase space region of the Pi and P3 
invariant mass close to the 1Z mass. Blob in the left diagram presents the non-resonant 
weak decay mode (see Fig. 1). 



APPENDIX 

Following the approach of the total amplitude contributing to the partial decay 
width T p can be written as a sum of the resonant and nonresonant contributions as defined 
in Eq. (21) in the following form: 

M = M nr + M r = Te~ h + P + R, (29) 

where T is the nonresonant tree contribution, P the nonresonant penguin contribution 
and R the resonant contribution to the amplitude. The partial width asymmetry defined 
in Eq. (22) is proportional to: 

A p oc sin 7 5s(T(P* + R*)) , (30) 

where 5s (A) stands for the imaginary part of A (similarly dt(A) stands for the real part 
of A). If we neglect the small imaginary part of the penguin Wilson coefficients, T and 
P will have the same strong phase. This implies that the only contribution to the partial 
decay asymmetry will come from the interference of the tree nonresonant and the resonant 
amplitude. One can write: 

A p oc $f(T)Vt(R) - 5s{R)^{T) . (31) 

The imaginary part of T is given by the absorptive part of the left diagram on the Fig. 
3. Using Cutkosky's rules, its contribution can be written as: 

5s(T)M(R) = J &(T)&(R)$l(S)d$ , (32) 

where the integration is taken over the P2P3 phase space. Here S denotes the strong 
re-scattering amplitude of P2P3 trough the resonance 1Z visualized in Fig. 4. Similarly, 
the imaginary part of R is given by the absorptive part of the right diagram on the Fig. 
3, where now the sum of all possible intermediate states into which 1Z decays should be 
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Figure 4: Re-scattering of the P2 and P3 states trough the resonance 1Z. 

taken into account. One can separate this contribution into the part with P 2 and P3 as 
an intermediate state ($s(R)p 23 ) and the part with all other intermediate states ($s(R)'). 
Again with the use of Cutkosky's rules, one obtains: 

5s(Rp 2 3 )^(T = J $t(T)&(R)$l(S)d$ . (33) 

The right hand sides of (J33)) and (J32j) are equal and therefore this two contributions to 
(13 lj) cancel among themselves and we have: 

A p ~ 5s(R)'$t(T) ~ F R (1 - Br(R - P 2 P 3 )3*(T) . (34) 

This cancellation is obviously a result of the unitarity and it maintains the equality of the 
total decay widths for the meson and the anti- meson as required by CPT theorem [HP]. 
That was already noticed by jHj and [T5] . where the more general proof is presented. 
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